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We present the experimental results on a series of high molecular weight, entangled polystyrene so- 
lutions subjected to a "mixed" shear and elongational flow-type generated in a co-rotating two-roll 
mill. For steady-flows as well as for start-up of flows from rest, we used two different optical methods 
to study the dynamics of entangled polymeric fluids: two-color flow birefringence (TCFB) and dy- 
namic light scattering (DLS). Using the TCFB method, we measured the birefringence, An, and the 
orientation angle, Xi of the optic axes in the solution and thereby the generalized viscosity function, 
rj (with the use of stress-optical relations). The DLS method was applied to measure the velocity- 
gradient, 7 , and the flow- type parameter A for the polymer solutions under flow conditions identical 
to the TCFB measurements. For low deformation rates the symmetry of the flow-field was reduced 
with the use of polymeric fluids compared to that seen with a Newtonian fluid. A molecular con- 
stitutive Doi-Edwards-Marrucci-Grizzuti (DEMG) model which includes polymer chain-stretching 
effects, has been used to numerically simulate the predictions for An, x, and 77 under steady flow 
conditions for the polystyrene fluids with 7 and A, measured via DLS, as inputs to the model. A 
detailed comparison of the DEMG model predictions with the experimental results shows that the 
model works qualitatively as well as quantitatively for the low and intermediate deformation rates, 
but fails at high rates of deformation by predicting a stronger chain-stretching than observed ex- 
perimentally. The effect of polymer molecular-weight and number of entanglements per chain are 
highlighted. 



I. INTRODUCTION 

The majority of prior experimental studies of the dy- 
namics of polymeric fluids were based on the simple shear 
flow. The main advantage is that the complete kinemat- 
ics of these homogeneous flow-fields is known a priory, 
and experimental measurements of velocity-gradients are 
not necessary. On the other hand, the majority of the 
polymer processing applications employ flows that are 
inhomogeneous and often involve a mixture of simple 
shear and pure extension. This have provided one moti- 
vation for recent experimental as well as theoretical stud- 
ies of the behavior of the polymeric fluids in "mixed" or 
purely extensional flows (also called the "strong" flows 
which in two-dimensions would require the magni- 
tude of the strain-rate to exceed the vorticity). Gener- 
ally, such strong flows appear only as a local part of some 
otherwise globally inhomogeneous and "weak" [Q flows. 
Examples of such flows include the stagnation zone in 
a cross-slot or cross-jet device; the converging zone of a 
contracting channel; and the stagnation region in two- 
or four-roll mills. The velocity-fleld in these local re- 
gion for any polymeric fluid will be very different from 
that of a Newtonian fluid of similar density and viscos- 
ity, because the flows are globally inhomogeneous and 
the fluid's viscoelastic nature. It is therefore insufficient 
to measure only the birefringence (or stress in a rheolog- 
ical experiment), without simultaneously measuring the 
velocity-gradient field by some technique. 

Different indirect techniques, e.g.. Particle Imaging Ve- 
locimetry and Laser Doppler Velocimetry pi, hot wire 



anemometers ||^, forced Rayleigh scattering [J], holo- 
graphic grating methods [^, etc. have been employed 
in the past for the measurement of velocity-gradients 
in time varying flows, but each exhibits some charac- 
teristic drawbacks. Provided that the flow is laminar 
and nonchaotic, the dynamic light scattering (DLS) tech- 
nique 1^,0 is one extremely valuable method for the di- 
rect measurement of velocity-gradients in strong flows. 
The high degree of spatial resolution of the laser light as 
well as the very high temporal resolution of the present 
day correlators has made it possible to use DLS for pre- 
cise determination of velocity-gradients for both steady 
and time-dependent flows even in very small-scale flow 
systems with velocity-gradients that change relatively in 
time. Further, unlike other methods mentioned above, 
DLS is nonintrusive in nature. 

The most frequently used indirect rheo-optical probe of 
polymer conformation in two-dimensional flows has been 
the flow birefringence technique Since it provides an 
average measure across the flow in third direction, it is 
only approximately local in that sense. Early studies us- 
ing traditional flow birefringence technique were primar- 
ily restricted to steady flows. On the other hand, recent 
studies [|-13| for both dilute and concentrated polymer 
solutions in time-dependent flows have used two-color 
flow birefringence (TCFB) iQ and/or phase- modulated 
birefringence techniques to simultaneously measure 
the retardance (and hence the birefringence, An) and the 
orientation angle, x, of the principle axis of the refractive- 
index tensor with respect to the axes fixed on the flow- 
cell, in a single experiment. 
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The theoretical understanding of the dynamics of en- 
tangled polymer solution and melts is far from complete. 
Using de Gennes' [|l6| original idea of "reptation", Doi 
and Edwards (DE) proposed a tube model which 
views the polymer chain as a sequence of segments con- 
fined in a tube-like region. The tube is defined by the 
topological constraints on the lateral motion and orien- 
tation imposed by neighboring chains. The segmental 
stretching can relax on a timescale similar to that for an 
unconstrained chain (dilute systems), i.e., the Rouse time 
T/j, and the segmental orientations can only relax via rep- 
tational dynamics i.e., through longitudinal diffusion to 
escape the tube which requires a much longer timescale 
Td, called the reptation or disengagement time. These 
two time scales are related via N^, the number of entan- 
glement points per chain: Td = SNeTn. The DE model 
is expected to be vahd for highly entangled {Ne >> 1) 
polymeric samples. Thus the original DE model neglects 
the chain-stretching effects by assuming that the "snap- 
back" of the stretched chain to its equilibrium length is 
instantaneous (i.e., t/j << Td) and develops constitutive 
equations which incorporate segmental orientational dy- 
namics only. This simplified model has been quite suc- 
cessful in many rheometric flows such as steady simple 
shear and oscillatory shear, where chain-stretching is ex- 
pected to be insignificant. On the other hand, serious 
inherent limitations of DE theory have long been recog- 
nized for the case of strong, extension-dominated flows 
P- p^ , p^ and also in transient shear flows, where chain- 
stretching is important. Over the past decade, this model 
has seen several improvements to include chain-stretching 
non-uniform 



in the uniform |11 
extensible 
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21| and finitely- 



p2| form. The predictions of the resulting so- 
called Doi-Edwards-Marrucci-Grizzuti (DEMG) model 
has been studied in detail in a recent set of papers by 
Mead et al. [||,|||. 

The aim of the present paper is to understand the dy- 
namics of entangled, high molecular weight polymeric flu- 
ids subjected to the strong, extension-dominated flows 
(i.e., for strong two-dimensional flows in the sense of As- 
tarita To this end, we present two-color flow bire- 
fringence and dynamic light scattering results for a se- 
ries of entangled polystyrene solutions in steady, near- 
homogeneous, planar flows created by a co-rotating two- 
roll mill. A detailed comparison of the experimental 
results for birefringence, orientation angle and general- 
ized extensional viscosity (defined below, and obtained 
via stress-optical relation) for the three entangled sam- 
ples to the predictions from the numerical simulation of 
the DEMG model ||2^ is carried out using the measured 
flow data as input to the model. The effect of molecu- 
lar weight and number of entanglements per chain of the 
polymers on the results are demonstrated. 

In sec. II we provide the necessary details related to 
our experiments, namely, the samples we used (sec. II. A); 
the linear viscoelastic measurements (sec. II. B); and the 



experimental apparatus (sec. II. C). The two-color flow 
birefringence and the dynamic light scattering techniques 
and apparatus involved are described in brief along with 
the two-roll mill flow-cell in sec. II. C. Section III deals 
with the results of the steady-state flow experiments 
(sec. III. A) using the DLS method and comparison of 
the steady flow TCFB results with the DEMG model 
(sec. III.B). Finally, sec. IV contains a summary of our 
findings and conclusion. 



II. EXPERIMENTAL DETAILS 
A. Materials 

The Newtonian fluid sample used was a suspension 
of spherical polystyrene particles (polyballs) of diame- 
ter 0.11/im (Polysciences Inc., U.S.A.) in glycerol (Sigma 
Chemical Co., U.S.A., ACS Reagent) at a concentration 
of 150 ppm. As shown in Table |, three non-Newtonian 
fluids, named PS81, PS82 and PS2, were made us- 
ing two standard, fairly monodisperse polystyrene sam- 
ples (Tosoh Co., Japan), namely = 8.42 x 10^, 
M^jMn = 1.17, Lot. No. TS-31 for one and M^, = 
2.89 X 10^, M^/Mn = 1.09, Lot. No. TS-10 for the 
other. Here and M„ specify the weight-average and 
the number-average molecular weights respectively. The 
solvents were prepared by adding polystyrene oligomer 
with Mjjj = 6000 or 2500 and a broad molecular weight 
distribution to toluene (Aldrich Chemical Co., Inc., 
U.S.A., ACS spectrophotometric grade) at mixing ratios 
by weight (toluene: 2500 PS = 43:57 for PS81 and PS82, 
and 48:52 for PS2), and allowing the mixture to dissolve 
at room temperature (= 20° C) for a few days with oc- 
casional slow mixing using a magnetic stirrer. Measured 
amounts of high molecular weight polystyrene samples 
were then added to these oligomer solvents to achieve 
the desired concentrations, c, shown in Table ||. The 
complete solution was then thoroughly mixed in a glass 
bottle for several hours by continuously rolling the bottle 
on its side in a RoUacell (New Brunswick Scientific Co., 
Inc., U.S.A.). In order to reduce the background scatter- 
ing in the light scattering experiments, the samples were 
filtered using Milipore filters with 5/im pore size for the 
polymer solutions and 0.45/im pore size for the polyball 
suspension. 

The flow-cell was pre-cleaned with acetone (Fisher Sci- 
entific, U.S.A., ACS spectranalyzed) , oven-dried at 50°C 
and brought back to the room temperature everytime be- 
fore the fluid sample was transferred into the cell via an 
air-tight siphon system. To minimize the loss of toluene 
in the process of fluid transfer, the use of an inactive gas 
in the siphon system was necessary. For this purpose, we 
used compressed nitrogen gas to maintain a minimum 
pressure differential in the siphon to fill up the cell vol- 
ume (~ 100 mL) in a couple of hours. 
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FIG. 1. The modulus of complex dynamic viscosities, 
(in P), as a function of the radial frequency, lo (in rad/s), for 
the three polystyrene samples, measured at 20° C. 
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FIG. 2. The dependence of the elastic modulus, G', and the 
viscous modulus, G", (both in the units of dyn/cm"^) on angu- 
lar frequency, ui (in rad/s), for the three polystyrene samples 
at 20° C. 



B. Linear viscoelastic measurements 

The linear viscoelastic shear flow properties of the 
three polymer solutions were measured, at a temperature 
of 20° C, using a Rheometrics Mechanical Spectrometer 
(RMS-800) in a cone-and-plate geometry with a cone di- 
ameter of 40 mm and a gap angle of 4° . The measured 



modulii of complex dynamic viscosities \ri*\ of the three 
liquids are plotted against the angular frequency to in Fig. 
|l|. We note that the viscosity of PS2 is about an order 
of magnitude higher than the other two samples; it is a 
much less mobile liquid. The sample PS81 is a couple 
of times more viscous than PS82. The typical nature of 
the three curves seen in Fig. |^ is very similar, namely, 
1 77* I is almost constant at low frequency a;<3xl0~^ 
rad s~^ and beyond this uj it shows a "shear" thinning 
behavior, following an empirical power law \r]*\ — Klo"', 
over the entire experimental frequency range, with the 
correlation coefficient being about 0.99 or higher in 
all three cases. The power law constant K and the expo- 
nent n, obtained from a linear regression of the data are 
given in Table |. Wc note that although the values of K 
are significantly different in these samples, the slope of 
the complex dynamic viscosity versus angular frequency 
curves in a double-logarithmic plot in Fig. Q is about 
same for two solutions with Ne ^13 and larger than the 
slope with Ne ^ 7. 

From the low frequency Newtonian plateau of these 
curves, the zero "shear-rate" viscosity rjo is extracted for 
these samples and is shown in Table ^. The number of 
entanglement points per chain, Ng, can be defined by 
analogy with the classical theory of rubber elasticity. 
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where c is the polymer concentration in g cm~^. Me = 
M^/Ne is the molecular weight between entanglements, 
R is the gas constant and T is the absolute temperature. 
We have used the approximate expression for the plateau 
modulus. 



G% = 3.44 X 10^ X c^-^ 



dyn/cm (c < 0.1 g cm '^), 

(2) 



obtained experimentally by Osaki et al. [£5| for 
polystyrene solutions in a good solvent. The values ob- 
tained using this equation agrees reasonably well with 
measured values. Thus Eqns. (0) and (0) leads to 



Ne = 1.41 X 10""* X M,„ X c} 



(3) 



for our experiments and the calculated values of A^e for 
the polystyrene samples are shown in Table |. 

The dynamic modulii of PS81, PS82 and PS2 are mea- 
sured on the same rheometer RMS-800 with the same 
pair of cone and plate as for the data in Fig. |^ and are 
shown in Fig. ||. Both, the elastic modulus (G") and the 
viscous modulus (G") show a strong frequency depen- 
dence for all three cases. The values of the modulii G' 
and G" for PS81 is higher than that of PS82, both hav- 
ing an order of magnitude lower values than in the case 
of PS2. Each of these pair of complex modulii shows a 
pattern of G" dominating over G' in the low frequency 
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"terminal zone" and the values of G' taking over that of 
G" in the higher frequency "plateau region", typical of 
polymeric fluids. In the low frequency flow region, the 
storage modulus G' approaches a quadratic dependence 
on uj and the loss modulus G" shows a linear depen- 
dence, again a typical behavior of polymers with a nar- 
row molecular weight distribution. At higher frequencies 
G' and G" tend to converge, which may be due to the 
formation of transient entanglement network The 
value of oj at which the storage modulus falls below the 
loss modulus is highest in PS2 and lowest in PS81. As- 
suming that this crossing point frequency corresponds to 
the inverse of the disengagement time r^, we have esti- 
mated the values of to be 67.50 s, 54.18 s and 16.80 
s respectively for PS81, PS82 and PS2. Using the pre- 
diction tr = Td/{iNf.) from reptation theory ||l^, we 
have calculated the Rouse time, t/j, for each sample as 
reported in Table ||. To check these values, we can use 
the correlation of Menezes and Graessley 
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where for polystyrene we use {Mc)q — 33, 000, 

a = 3.4, b — 1.3, density p — 1.07 g cm~^ and volume 
fraction v = pc. The calculated Rouse times are 1.95 s, 
3.84 s and 0.87 s for PS81, PS82 and PS2 respectively. 
It is unclear to us why the values calculated from Eqn. 

are not in closer agreement with those obtained here 
experimentally, though they are in the same ballpark. 
In what follows, we use the experimentally determined 
values listed in Table ||. 



C. Experimental set-up 

A schematic diagram of our experimental system is 
shown in Fig. ^. It is composed of two main parts: the 
TCFB optical arrangement that is used to measure the 
optical anisotropy of the fluid and the DLS set-up that 
is used to measure the flow parameters of the fluid. In 
the following, we will briefly touch upon some specific 
details of the apparatus relevant to the TCFB and DLS 
studies reported in this paper. Further details of the data 
analysis technique used for TCFB ||l^ and DLS have 
been described elsewhere. 




FIG. 4. Schematic representation of the light scattering 
configuration. 
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FIG. 3. Relative orientation of the optical system including 
the flow-cell and the definition of the coordinate systems. 



1. Two-color flow birefringence 

The TCFB technique [|lj] simultaneously measures 
both retardance and orientation x of the principal axis of 
the refractive- index tensor n relative to the axes [{x, y) in 
Figs. H andlJ] fixed in the flow device. The light source is 
an Argon-ion laser (Spectra Physics Model 2020) operat- 
ing at 300 mW, which emits two intense wavelengths at 
Xb = 4880 A(blue) and Ac = 5145 A( green). As shown 
in Fig. 1^, with the use of several optical elements, each 
of these beams (independently polarized at 45° relative 
to each other) are made to pass along an identical opti- 
cal path through the sample in a two-roll mill which is 
placed between crossed polarizers. Since the flow is non- 
homogeneous, to ensure the collinearity of the two beams 
(with identical optical properties, e.g., beamwidth, Gaus- 
sian beam-profile etc.) passing through the same element 
of fluid in the flow-cell is very crucial in this set-up. Each 
of the measured intensities (at the photodetectors) de- 
pends upon the current degree of optical anisotropy of 
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the sample in the plane of flow [{x,y) plane in Fig. |[, 
— i.e., the birefringence, An = ny — n± — and %, which 
provides a measure of the average orientation of Kuhn 
segments in the polymeric liquid. Further details of the 
optical system are described in Ref. and will thus 
not be repeated here. The measured overall error due to 
nonidealities of the optical components indicate that the 
maximum extinction ratios detectable in this set up with 
the flow device (loaded with polymer solution), for both 
colors, are typically 0(2 x 10~^). 
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FIG. 5. Schematic showing the optical layout for the 
two-color flow birefringence and the dynamic light scattering 
setup. M: mirror; F: interference filter; L: lense; P: polarizer; 
A: analyzer; D: photodiode; S: pinhole; 6: scattering angle. 
The subscripts B &l G specify the blue (continuous line) and 
the green (broken line) laser beams, respectively. The super- 
scripts o and e refer to the ordinary and extraordinary beams, 
respectively. 



2. The two-roll mill and the flow 

The two-roll mill is a flow device consisting of two 
cylinders driven simultaneously by a single DC stepping 
motor (Superior Electric, U.S.A., SLO-SYN motor type 
MO62-FD-09). The cylinders rotate at the same angu- 
lar velocity and in the same direction to generate a flow 
between the roller pair, that can be approximated at the 
central stagnation point by a linear (or homogeneous), 
planar velocity-field of the form 



V — Vv • r. 



(5) 



Here, r is the position vector defined in the flow plane 
w.r.t. the {x,y) frame, shown in Figs. |^ and ^ The 
velocity-gradient tensor Vv in this coordinate system is 
given by 



Vv — 7 



e 1 
A -e 



< A < 1. 



(6) 



If we assume further that the flow is symmetric about 
the central {x, z) plane between the rollers, as in the case 
for a Newtonian fluid at zero Reynolds number, then 







(symmetric flow). 



(7) 



In that case, Eqn. (^ can be completely characterized 
by two scaler parameters, namely, the magnitude of the 
velocity-gradient 7 = |Vv| (also called the "shear-rate") 
and the flow- type parameter A. The parameter A is a 
measure of the ratio of strain-rate to vorticity, and is 
deflned as, 



1^1 1 
IIOI 



A 



1- A' 



(8) 



Thus, A = corresponds to a simple shear flow, A = 1 
to a purely extensional (also called hyperbolic) flow, and 
intermediate values of A represent a "mixed" shear and 
elongational flow-type (also called a "strong" or exten- 
sion dominated flow , since the strain-rate exceeds the 
vorticity). In sharp contrast to other extensional flow 
devices, polymer molecules at the stagnation region of 
a two-roll mill have very long residence time and hence 
are subjected to very large total strains. Consequently, 
they can, in principle, achieve the maximum change in 
conformation that is consistent with a particular strain- 
rate For a Newtonian fluid, at very low Reynolds 
numbers, where the symmetric flow assumption is valid, 
the creeping flow solutions Q for a two-roll mill in an 
unbounded fluid approximately relate the values of the 
corresponding velocity-gradient 7jv and the flow-type pa- 
rameter Aat in the stagnation region with the geometry 
of the cell (i.e., the roller radius R and the gapwidth 2/i), 
via 
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with A — 



Au) 



K 



f 



Kf + sinh cosh Kf ' 



(9) 



where oj is the angular velocity of the rollers and the 
parameter Kf is given by 



1 



h 
R 



cosh K 



(10) 



The symmetry axes of the flow- field, the principal optical 
axes of the solution, as well as the relative orientation of 
the blue polarizer and analyzer for the two-roll mill set 
is shown in Fig. ^. It may be noted that the Newtonian 
value of the flow-type parameter, Aat, is also related to 
the acute angle of crossing, 2(p, of the streamlines passing 
through the stagnation point, according to 



A 



N 



tan tp. 



(11) 



5 



The present two-roher can, in general, be used [|To[ with 
a pair of rollers chosen from a set of eight such pairs of 
different diameters covering < A < 0.25. For this study, 
the dimensions of the rollers used are specified in Table ||, 
which corresponds to A = Xn = 0.1501 for a Newtonian 
fluid. The subscripts "th" and "exp" in Table || refer 
to the theoretical and experimental values, respectively. 
Eqns. (||) and (10) allow us to estimate the theoretical 
value of (ji^/uj)^-^ for this set of rollers, as given in Table 

The stepping motor is interfaced to a computer 
(Hewlett Packard, model 9133 for TCFB or PC-486 for 
DLS) via a GPIB (National Instruments, U.S.A., model 
NI-488.2) switch board as shown in Fig. |5[ The motor 
has a fast response time of C(10 ms) and a maximum 
acceleration of 100 000 steps s~^. Gears with reduction 
ratios 5 : 1 or 20 : 1 are used between the motor and the 
flow device to cover a wide range of velocity-gradients. 
The accuracy involved in repositioning of the flow cell 
with respect to the incident beam, even when it is dis- 
mounted for the purpose of replacing the solution, is al- 
ways better than 0.0025 cm. in x, y and z directions (as 
defined by three translation stages used to control the 
respective movements) and 0.001 rad in the azimuthal 
orientation of the flow-cell (defined by the micrometer 
used to control the orientation) [Fig. The flow-cell was 
thermostated within ±0.2°C via a temperature-regulated 
waterbath circulator. 



3. Dynamic 



scattering 



The optical setup for the dynamic light scattering ex- 
periment was built around the TCFB apparatus, as can 
be seen from Fig. ^ by mounting the necessary opti- 
cal accessories on a triangular optical rail placed on the 
rotating arm of a massive goniometer that defines the 
scattering angle, 9, and hence the scattering vector q 
[q = \q\ = ^^sm{0/2) as shown in Fig. ^, where n is 
the refractive- index of the solution] . The green beam of 
the laser is blocked using a beam stop near Pq, in Fig. |5[ 
The incident blue beam scattered by the sample is polar- 
ized and coUimated to project an image of the scattering 
volume at the photomultiplier tube (PMT) [Hamamatsu, 
model R647-04]. A blue line filter, Fb, obstructs any spu- 
rious light from reaching the detector. The preamplified 
and discriminated PMT pulses are fed to a 72-channel 
correlator (Brookhaven Instruments, U.S.A., model BI- 
2030) to construct time-resolved intensity autocorrela- 
tion functions. The commercial correlator software has 
been modified and used in PC-486 computer to control 
the flow experiment, as well as to analyze the correlation 
functions 0. 

Provided that the seed particles are isotropic scat- 
terers, and the time scale associated with the velocity- 
gradient [t-y {qjL)~^, where L is the laser beamwidth] 



is much smaller than the time scale for diffusive motion 
[to ^ iDq'^)~'^], the homodyne intensity correlation func- 



tion for a general linear flow in the beam coordinates 
ix',y',z') [Fig. I is given by 



F2{q,t)^P' 



(fr'I{f') exp{~iq' ■ Vv' 



r't} 



(12) 



Here, /?' is the coherence factor defined by the optical 
geometry . The velocity-gradient tensor Vv of Eqn. 
(^) can be expressed in the beam coordinates as 



Vv' = Q • Vv • where Q 



cos ( 
sinf 



— sm ( 
cosdi 



(13) 



In Fig. |[ the beam coordinates are chosen such that the 
scattering vector (f is orthogonal to the y' axis, i.e.. 



g'- Vv'-r' = <7,/(0)x' + g,g((/))y', 



(14) 



where qx ~ qcos{9/2), f{(f>) = cos0(ecos0 — jXaincf)) — 
sin0(7cos0-|-esin(/)) and (;((/)) — sin(/)(ecos0— 7Asin(/))H- 
cos (/)(7 cos (/) -I- e sin (/)) . 

Using a Gaussian intensity profile, Eqns. ( |l3|) and ([ij) 
in Eqn. (02), we get 



i^2(g,t) =/3exp^ --iVi'cos^ I - ) h{<i>) \ , (15) 



where /3 oc csc^ Q and 



M0) = /'(0)+5'('/')- 



(16) 



The values of /i(0) are evaluated in Table |T| for three dif- 
ferent azimuthal orientations of the two-roll mill. Thus, 
with the assumption of a symmetric flow [Eqn. Q , 7 can 
be obtained from the decay rate /i(0) of the measured cor- 
relation function [Eqn. (|l5|)] at one orientation, = 0, 
of the flow-cell and A from an additional measurement at 



a second orientation, — 90° (see. Table HI). The cor- 
relation functions were accumulated only after the flow 
had attained its steady value for the corresponding mo- 
tor speed. To improve the signal to noise ratio, each cor- 
relation function was obtained by averaging over many 
repetitions of the experiment, ranging from a minimum 
of 200 repetitions at the high roller speeds, up to a max- 
imum of 800 repetitions at the low roller speeds. Using a 
simulated annealing Monte Carlo fitting procedure [p9| , 
these average correlation functions were then fitted to 
Eqn. (^). In order to verify that the experimentally ob- 
tained correlation functions are very close to exponential 
in nature so that the above procedure followed to extract 
the velocity-gradient components from the decay time is 
indeed justified, we required a correlation coefficient i?^, 
specifying the quality of the fit, exceeding 0.99. 

Contrary to the basic assumption of optical isotropy 
of scatterers used to derive Eqn. (|l^), Wang et al |^] had 
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shown that this equation could also be used for the scat- 
tering from (intrinsically anisotropic) polymer molecules. 
In that case, the pre-exponential factor (3 is shown to be 
directly related to the components of the intrinsic molec- 
ular polarizability tensor. 
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Position X (mm) 

FIG. 6. The velocity-gradient profile along the central 
crossing line of the rollers with A = 0.1501 and = 0° for 
the Newtonian fluid at w = 1.07 rad/s (triangles and circles) 
and uj = 1.57 rad/s (diamonds). The stagnation point was 
defined as the zero position. The hollow symbols denote the 
measured data and the filled symbols denote the data ob- 
tained by symmetry. 



III. RESULTS 
A. DLS steady flow results 

1. Flow characterization 

The DLS technique provides us the opportunity to map 
the flow-field in a large region between the rollers. In or- 
der to look for changes in the steady flow-field due to 
viscoelasticity in the presence of polymers compared to 
that seen with a Newtonian fluid, we first characterize 
the flow using a Newtonian fluid. This provides a com- 
parison with predictions from the creeping flow theory 
for an unbounded two-roll mill. The measured velocity- 
gradient 7 along the x axis through the stagnation point 
of the two-roll mill orientated at (p — 0° for an angular 
velocity uj = 1.57 rad/s is plotted along with two sets of 
data obtained by Wang et al Q, in terms of 7/ci; in Fig.|6[ 
The filled symbols denote the points where the geometric 
construction of the flow-cell did not allow measurements, 
and those points are therefore obtained by the symmet- 
ric reflection of the hollow measured points at negative 



X values. As can be clearly seen from this figure, differ- 
ent sets of experimental results obtained at widely sep- 
arated times with different roller rotation rates overlap 
perfectly. In a small region surrounding the stagnation 
point the velocity-gradient were approximately constant, 
7/cj = 2.8, certifying that the flow-field in this region 
is nearly homogeneous. Also, this value compares very 
well with the theoretically expected value of 7Ar/w = 2.7 
(Table. |l|). From here onwards, we shall use 7jv and 
in the text to specify the values of 7 and A [Eqns. (|^) 
and ( p^ j at the stagnation point of a two-roll mill filled 
with a Newtonian fluid, and the imposed motor speed 
will be measured in terms of the corresponding Newto- 
nian strain-rate 'JnV^n ■ 
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FIG. 7. Linear dependence of the magnitude of the veloc- 
ity-gradient components, 7 and 7A (both in the units of s~^ 
and measured at </> = 0° and (j) = 90°, respectively), on the 
Newtonian rates of deformation, Jn^]^'^ (in s), for a partic- 
ular set of the rollers, Xm ~ 0.1501, in steady flow at the 
stagnation region for the Newtonian fluid. The inset shows 
the strain-rate dependence of the flow-type parameter. A, ob- 
tained by dividing one set of data by the other. The straight 
lines represent the creeping flow solutions. 

In Fig. 1^, we show the theoretical predictions for 
and t-atAat [Eqn. (^), (10) and Table |l| versus jNVXj^ 
in the form of straight lines. The experiments were car- 
ried out for 14 different Newtonian strain-rates, and for 
both parallel and perpendicular orientations of the sam- 
ple cell. To extract the experimental values of 7 and 7A 
from the measured decay rates of the correlation func- 
tion of Eqn. (|l^) at (j) — 0° and 90° respectively, we 
required a value for the beamwidth L. With an initial 
guess of L = 29 /Ltm, that was reported for a previous set 
of experiments |^ from this laboratory, we used a sim- 
ulated annealing Monte Carlo fitting technique 123] to 
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find via the fit shown in Fig. ^ that L = 32 /^m yields 
the best match of the theoretical straight lines for both 7 
and 7A over the entire range of experimental strain-rates. 
Compared to the previous experiments , we have per- 
formed the present experiments including much higher 
motor speeds. It is clearly seen in the figure that the 
corresponding data agrees well with the creeping flow so- 
lution. The inset shows that the extracted value of the 
flow-type parameter A also maintains a constant value as 
expected. 
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FIG. 8. The dependence of the velocity-gradient, 7 (in 
s~^), on the Weissenberg number WiR = 'yX^^^ra (based on 
the Rouse time, th and the measured 7 & A), for Ajv = 0.1501, 
in steady flow at the stagnation region for the three entangled 
polystyrene fluids. In the inset, the same data are plotted ver- 
sus the Newtonian rates of deformation, 7jvA]^^ (in s~^), and 
are compared with the creeping flow solution for 'yN (straight 
line) . 

In a similar manner, the flow parameters are mea- 
sured for the three polymeric samples, PS81, PS82, and 
PS2, subjected to steady-state flow conditions. The re- 
sults are plotted in Figs. ||, §, and |l^ versus the non- 
dimensionalized deformation rate. Win — 'yy/Xrn , called 
the Weissenberg number (based on the Rouse time, r^, 
and measured values of 7 and A). The Weissenberg num- 
ber is defined as the ratio of a characteristic relaxation 
time for the polymeric fluid to a characteristic time for 
deformation. It specifies the ability of a flow to generate 
departures of the polymer configuration from its static 
equilibrium value. The same data are also plotted against 
JnVXn , in the inset of these figures, for a comparison 
with the Newtonian values shown by the straight lines. 
As expected, there is a substantial reduction in both 7 
and 7A from the corresponding Newtonian value. The 
deviation from the Newtonian values increases with the 



increasing roller rotation rate. Somewhat unexpectedly, 
however, both 7 and 7A still appear to increase approxi- 
mately linearly with the Newtonian strain-rate. It is also 
noteworthy that this increase is very similar for the two 
samples, PS81 and PS2, which have the similar number 
of entanglements (TVg ^ 13) per chain (see. Table ^. On 
the other hand, PS82, with ^ 7 shows a different slope 
for its almost linear increment with increasing 'yN\^N ■ 
When plotted against the Weissenberg number (Figs. || 
and 1^) , the qualitative feature of the curves for 7 and 7A 
appear to be quite different than that against JnVXn 
: they become much more linear and owing to an order 
of magnitude smaller r^j, the increment for PS2 with the 
increasing rate of deformation is much steeper than the 
other two samples. Also, 7 versus Win data for the two 
samples with similar molecular weight but different Ne 
are similar up to Win ~ 3. The quantitative similar- 
ity in deformation rate dependence for the samples with 
similar is more apparent for the flow-type parameter 
A, as can be seen in Fig. |o[ The qualitative nature of 
this dependence for all three samples appears to be sim- 
ilar, namely, it reaches an almost constant value (~ 0.1 
for the samples with 13 and ~ 0.05 for PS82) at 

high deformation rates Win > 0.4 (7ivV^Af > 1) and 
at low rates of deformation exceeds the Newtonian val- 
ues. There is an intermediate transition regime between 
these these two different behaviors. The parameter A 
attaining values higher than the Newtonian value at the 
lowest Win studied, seems surprising since in the limit of 
extremely low deformation rate (Win — s- 0), one should 
expect even a Non-Newtonian fluid to exhibit Newtonian 
behavior. We are still away from this limit in that the 
Weissenberg number corresponding to the lowest motor 
speeds for the steady flow experiments are ~ 0.05 (PS2), 
~ 0.1 (PS81) and - 0.2 (PS82). As mentioned earlier, 
each data point in 7 and 7A curves are obtained by aver- 
aging over several repeated experiments and then the A 
values are obtained by point to point division of the two. 
Following a similar procedure, we have calculated the 
error-bars on the A values from the standard deviations 
of the repeated experiments for 7 and 7A and plotted in 
Fig. As expected, the error-bars increase for lower 
Win, because of the division of smaller values of the cor- 
responding 7A by 7 . Even within the limits of these 
error-bars we can clearly see that A initially exceeds Ajv ■ 
We have confirmed that this finding is not an experimen- 
tal artifact by two other experimental means employed 
in this study, namely, the measurements of birefringence 
and flow parameters in transient flow conditions, and is 
described in section III. A. 3 of this paper. 
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FIG. 9. The dependence of the component of veloc- 
ity-gradient measured at = 90°, i.e., 7A (in s^^), on the 
Weissenberg number Win = 'yX^^^TR, for Ajv = 0.1501, in 
steady flow at the stagnation region for the three entangled 
polystyrene fluids. In the inset, the same data are plotted 
versus the Newtonian rates of deformation, ijvAj^^ (in s~^), 
and are compared with the creeping flow solution for 7jvAjv 
shown by the straight line). 
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FIG. 10. The flow-type parameter A versus Win for 
Ajv = 0.1501, in steady flow for the three entangled 
polystyrene fluids, extracted by dividing the data in Fig. 8 by 
the corresponding data in Fig. 7. The inset shows the same 
data, plotted against 7jvAJ/^ (in s~^) and compared with the 
creeping flow solution of Xn = 0.1501 (straight line). 



WiR 

FIG. 11. The error-bars represent the standard deviation 

on the flow-type parameter, A, measured at several Weis- 
senberg numbers, Win, for the steady-state flow of the en- 
tangled polymeric solutions. 




Calculated 7(1+^^)'* (s ') 
FIG. 12. The connected symbols illustrate the dependence 
of the velocity gradient component measured at (p = 45°, i.e., 
7(1 + A^)^''^, against its expected value, calculated using 7 
and A values measured at = 0° and = 90°, respectively, 
for the three entangled polystyrene fluids in steady flow at the 
stagnation region. In a similax fashion, the experimental and 
the theoretical (creeping flow) data for the Newtonian fluid 
are shown by the hollow squares and the solid line, respec- 
tively. 



9 



2. Flow symmetry 

As we have mentioned before, the assumption that the 
flow in the two-roll mill is symmetric about the plane 
passing exactly midway between the co-rotating rollers 
{x, z) plane in Figs. ^ and ^) is valid for the creeping 
flow of a Newtonian fluid. For a viscoelastic fluid, one 
might expect the symmetry of the flow to change com- 
pared to the Newtonian case. Using the measured val- 
ues of 7 and 7 A (Figs. || & ||) at several different steady 
Newtonian strain-rates, we have calculated the expected 
values of 7-^/ (1 + A^) . Then, by repeating exactly same 
steady flow conditions that were used to measure 7 and 
A, we have extracted the velocity-gradient component 
7"^ (1 -I- A^) from the decay rates of the correlations func- 
tions measured directly with the two-roll mill oriented at 
(j) = 45° (see, Table III). The result is plotted versus 
the corresponding calculated values at the same Newto- 
nian strain-rates with the connected hollow symbols in 
Fig. |2| for the three viscoelastic samples. The hollow 
squares represent the measured values of 7-^/ (1 + A^) for 
the Newtonian fluid a.t (j) — 45° . Both these and the the- 
oretical values, i.e., 7Ar-\/(l -I- A^) (shown by the solid 
line) at several strain-rates are plotted in Fig. |l^ against 
the expected values of the same parameter that are cal- 
culated using the measured 7 and 7A data from Fig. |^ 
If the symmetry of the flow is maintained, then Eqn. 
( p^ ) should be strictly valid for all orientations, (f>, of the 
flow-cell. This means that each curve in Fig. |l2| should be 
linear. As can be clearly seen from the figure, this is best 
followed for the Newtonian fluid indicating the fact that 
the flow is symmetric in this expected. On the 

other hand, the data for the entangled fluids indicate that 
the flow symmetry is maintained only approximately, i.e., 
within the experimental accuracy of 10%, for the first 
five or six data points (i.e., for "fNy^N < 0.5 s^^ ). The 
data for PS81 and PS2 (both with ~ 13) deviate more 
from the linearity with the increase of 'Jn^/^n , but for 
PS82 {Ne ^ 7) the experimental data are fairly close to a 
straight line over the entire range of strain-rates studied. 
At high rates of deformation, the measured 71/ (1 -I- A^) 
is higher than the expected value for PS81 and PS82, but 
is lower in the case of PS2. 

In order to quantify how flow symmetry changes with 
the deformation rate, we now proceed to evaluate the pa- 
rameter e of Eqn. (^). To do this, we can first extract the 
values for h{(j3 = 0°), /i((/) = 45°) and /i(0 = 90°) from 
the experimentally measured decay rates of the correla- 
tion function at the three aforementioned orientations of 
the two-roll mill. When these values are used in conjunc- 
tion with Table III , we get three equations involving the 
three unknowns 7 , A and e. We have used the Newton- 
Raphson's method to numerically solve these equa- 
tions for the case of a Newtonian fluid as well as for a 
representative case of the viscoelastic fluids of our study. 



namely, for PS2. The results are shown in Fig. |l^. Our 
experiments with a Newtonian fluid shows that e is zero, 
i.e., the flow is perfectly symmetric, for strain-rates upto 
about 2 s~^ . Then, as the deformation rate is increased 
the value of e increases almost linearly, though its value 
remains extremely small (^ 10"^) even at the highest 
strain-rates studied, justifying the fact that, for a New- 
tonian fluid the flow-symmetry is maintained. This also 
explains why the first seven data points {'JnVXn < 2 s^^ 
) for the Newtonian fluid in Fig. |l^ lies on the theoretical 
solid line and then slightly deviates from the line for the 
higher values of ^nVXn ■ For a non-Newtonian fluid, 
on the other hand, the flow-symmetry is maintained for 
only up to 'yN\^N < 0.5 s~^ . Our transient flow ex- 
periments with these polymer solutions are carried out 
in this range of deformation rates (which corresponds to 
WiR < 1) and will be reported in a subsequent paper 

isi. 




Newtonian fluid 
PS2 



FIG. 13. The parameter e versus the Newtonian 

1/2 

strain-rate, 7jvA^ , for a representative polystyrene fluid PS2 
and the Newtonian fluid. The inset shows the "zoomed up" 
result for the Newtonian fluid. 



3. Steady flow and flow-type parameter 

Let us now look into the effect of the polymer solu- 
tions on the steady-state values of the flow-type param- 
eter A at different angular velocities of the rollers. We 
have briefly referred to this point in Fig. |l^. Since, we 
have also performed extensive DLS and TCFB experi- 
ments on these samples subjected to transient flows , 
namely, the startup flows from rest, for several constant 
Weissenberg numbers {Win < 1), it would be worthwhile 
to check how the asymptotic values of A extracted from 
each of those different experimental techniques vary with 
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Win and also how they compare with the steady flow re- 
sults of Fig. From the transient DLS experiments, we 
can directly get the asymptotic values of A but for the 
TCFB experiments, we should follow an indirect method 
to calculate A. This can be done because for Win < 1 
the flow retains its symmetry, as we have shown before. 
The steady-state orientation of the principal optic axis at 
appreciable strain-rates is expected to approach the out- 
flow axis of the flow, or, in other words, X ~^ x' ~ [see, 
Fig. 1^. Hence, A may be calculated from the asymptotic 
orientation angle, x', of the flow birefringence data |31 
for the startup flows, via a relation similar to Eqn 
i.e., 
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A tan^ x'- 



(17) 




also show A values exceeding Aat at the lowest rates of de- 
formation studied, which again are higher than Wir = 0, 
where we should expect that A = Aat. In fact, for PS81, 
at the three lowest rates of deformations studied, the ex- 
tracted value of A from the transient experiments do show 
a trend to decrease after reaching a maximum. 
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FIG. 15. Same as in Fig. 14, but for the polystyrene sample 
PS82. 



B. The DEMG model comparison of the steady flow 
TCFB results 
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FIG. 14. The Weissenberg number {Win = ■jX^^^tr) de- 
pendence of the flow-type parameter. A, for the polystyrene 
sample PS81, extracted from three different means (as 
explained in the text) and plotted together along with 
Ajv = 0.1501. The Weissenberg number is calculated based 
on the measured values of 7, A and tr. 

We have presented the results in Figs. |lj, ^ and 
H for the samples PS81, PS82, and PS2 respectively 
and compared the results with the theoretical value of 
Aat = 0.1501. For each of these three polymeric samples, 
the overall qualitative match of the data as well as the 
quantitative match between the A values calculated us- 
ing three different experimental procedures, is quite sat- 
isfactory. The fair quantitative agreement between the 
asymptotic values of A obtained from the transient flow 
experiments and the steady flow data for A confirms that 
the startup experiments with the entangled samples have 
almost reached their steady values in the total predeter- 
mined evolution time of te — 30 s pl| ]. In agreement 
with the steady-state result, the transient experiments 



Probably one of the most used theoretical models to 
describe the dynamics of entangled polymers is the Doi- 
Edwards (DE) tube (reptation) model [|l^. The primary 
drawback of this model, as is clearly evident by its fail- 
ure to predict the polymer dynamics at high deformation 
rates jVX > I/tr, is the fact that the primitive chain 
is assumed to be inextensible. Subsequently, in an effort 
to improve its predictions, the Doi-Edwards-Marrucci- 
Grizzuti (DEMG) model [^,^ was developed to incor- 
porate chain-stretching into the original DE tube model. 
As noted earlier, the DEMG model contains two widely 
separated time scales, one for the relaxation of orienta- 
tion (r^) and a much shorter one for the relaxation of the 
stretch of polymer segments {tr). Thus, as the strain- 
rate jy/X is increased, segmental orientation takes place 
first when 7a/A ~ I/t^ or Wid = jVXTd = 0{1) and seg- 
mental stretching then "starts" at a much larger strain- 
rate when WiR = jV^TR — 0{1). Also, it is expected 
that segmental orientation will take place earlier in di- 
lute systems than in entangled systems. In agreement 
with the aforesaid expectations, a number of numerical 
studies [|2l[ have suggested that the DEMG model has 
more success over the DE model when compared to ex- 
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perimental results for startup of a simple shear flow at 
high deformation rates {Wiu > !)• Following this, in 
a recent numerical study ||2^,|2^, the DEMG model was 
developed further by "adding" a non-linear finitely ex- 
tensible spring force. 




FIG. 16. Same as in Fig. 15, but for the polystyrene sample 
PS2. 
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FIG. 17. The experimental steady-state birefringence data 
for the three entangled polystyrene solution plotted as func- 
tions of the Weissenberg number, Win = 'yX^^^TR, based on 
the measured magnitude of the principal eigenvalue, 'yA^''^, of 
the velocity-gradient tensor, Vv, and the Rouse time scale, 
tr. The birefringence is normalized with CG%, where G% is 
the plateau modulus, and C is the stress-optical coefficient. 
The data are plotted in the same scale as the corresponding 
DEMG model predictions in Fig. 18. 



z lo' 
O 

u 

"a 
< 



10-' 



DEMG model predictions 




Wi„ 



FIG. 18. The DEMG model predictions for the 
steady-state An/ {CG%), using measured values of both 7 and 
A, at several different strain-rates Win but the same set of 
rollers, Ajv — 0.1501, for the polystyrene samples. The model 
parameters rit, Ne, and tr axe given in Table I. 

For an entangled solution, birefringence measurements 
directly relate to changes in the conformation of the poly- 
mer molecules. As mentioned earlier, the complete char- 
acterization of the flow- induced anisotropy of a polymeric 
fluid requires determination of the birefringence, An, as 
well as the orientation angle, Xi of the principal axes of 
the refractive index tensor. In this section, we present the 
steady-state two-color flow birefringence results and com- 
pare them with the predictions of the DEMG model us- 
ing the flow parameters measured via DLS as input. We 
assume that the velocity- gradient tensor in ^ 75/im di- 
ameter birefringence measurement zone surrounding the 
stagnation point of the co-rotating two-roll mill can be 
approximated by Eqn. (H) and @). The values of the 
model parameters TVg, r^, and nt calculated for the three 
solutions are given in the Table |[ Here rit is the number 
of Kuhn statistical subunits in the chain based on the 
assumption of 10 monomers per subunit and is calcu- 
lated from the molecular weight of the polymer |3|] . The 
steady-state experimental results for the dimensionless 
birefringence. An/ (CG'f^), for the three entangled sam- 
ples are presented in Fig. |l^ as functions of Win. The 
birefringence. An, is scaled by the birefringence CG'pf 
(obtained by using the stress-optical law) that would be 
present at a stress level equal to the plateau modulus, 
G^. Here, C is the stress-optical coefficient which in the 
free-jointed chain model is defined as 



C 



27r 

45 



(n + 2)2 



(18) 



where a\\ and a± are the components of the intrinsic 
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molecular polarizability tensor a{t) along and perpen- 
dicular to the Kuhn statistical segment, and n is the 
bulk refractive-index of the medium. This is justified 
only up to a modest level (< 50%) of fractional chain- 
extension, so that non-Gaussian effects can be neglected. 
To scale the experimental birefringence, an approximate 
expression given in Eqn. (0) for the plateau modulus, G^^r, 
and the literature value B of C = 5 x 10~^° cm^/dyn 
for polystyrene solutions is used. The model predictions 
corresponding to the experimental curves in Fig. |l^ are 
displayed in Fig. |l^ (both on the same scale). Each in- 
dividual point on the predicted curves of Fig. ^ is the 
steady-state value obtained from separate runs of DEMG 
model numerical simulations using the measured values 
of 7 and A for each corresponding point on the experi- 
mental curves (Fig. |lj). The model clearly predicts |24) 
the existence of three flow regimes in the steady-state, 
namely, the linear viscoelastic regime seen at low rates 
of deformation Wid < 0(1), the non-linear viscoelastic 
regime seen at the intermediate deformation rates lying 
between Wid > 0(1) and Win < 0(1), and the highly 
non-linear viscoelastic regime seen at high deformation 
rates Win > 0{1). In each of these flow regimes, the ef- 
fect of three parameters, namely, nt, Ne, and A on the na- 
ture of deformation rate dependence of the birefringence, 
orientation angle, stretch and viscosity, as predicted by 
the DEMG model, has been discussed in detail in a recent 
publication by Mead et al. p4[ . 
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FIG. 19. The experimental steady-state orientation an- 
gle, X (in degrees), for the three entangled polystyrene so- 
lution plotted as a function of the Weissenberg number, 
WiR — 'yX^^^TR, using the measured values of 7, A and tr. 
The data are plotted in the same scale as the corresponding 
DEMG model predictions in Fig. 20. 
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FIG. 20. The DEMG model predictions for the 
steady-state orientation angle, x (in degrees), using measured 
values of both 7 and A, at several different strain-rates, WiR, 
for the three polystyrene fluids in a two-roll mill with a fixed 
set of rollers, Ajv = 0.1501. The model parameters nt, Ns, 
and Tr are given in Table I. 

In the linear viscoelastic regime, the slope of each ex- 
perimental curve on the double-logarithmic plots (Fig. 
p^ ) is indeed unity ||2^, similar to the predicted curves 
(Fig. |lj) . According to the predictions of Ref. , in this 
flow regime, the relative positions of the birefringence 
curves from top to bottom for these entangled samples is 
primarily determined by the values of Ne and A (curves 
with lower N^, and higher A will fall lower), while nt has 
no effect. As clearly seen from Figs. |l^ and the lower 
value of for PS82 keeps the An/ {CG'^) curve lower 
than those of the other two samples. The lower value of A 
for PS82 in this regime (Fig. |o|) would predict otherwise 
. The slight difference between the curves for PS2 and 
PS81 (both having ~ 13) seen both in the experimen- 
tal and predicted plots, is due to the slight difference in 
the A values between the two samples in this regime, as 
can be clearly seen from Fig. [l^. In this flow regime, al- 
though the initial low levels of segmental orientation, as 
predicted in Fig. |2^, is indeed seen for the experimental 
curves in Fig. |l^, the rate of decrease of x with increas- 
ing Wiu, as seen in the experiments for all samples, is 
much faster than predicted. Also, the orientation angles 
for PS81 at a fixed value of Win in this regime lies in 
between those for the other two liquids which again is in 
accord with the theoretical expectations (taking into 
account the effects of the corresponding values of nt , 
and measured A). 

When the second flow regime of non- linear viscoelastic- 
ity is approached, Wid = 0(1), the scaled birefringence 
curves are predicted to depart from linearity and begin 
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to approach a plateau, as seen in the DEMG results of 
Fig. |l8|. The width of the plateau is directly proportional 
to A'^e, and is therefore expected to be about half as wide 
for PS82 {Ne ~ 7) as for PS2 and PS81 (TVe - 13). The 
experimental curves in this regime do show non-linearity 
for all three samples, but the collapse of the curves and 
an existence of a narrow plateau is seen only for PS81 
and PS2 and not for PS82. The dynamics in this regime 
are controlled by strong segmental orientation of the tube 
without stretching. Referring to Fig. we see that the 
monotonic decrease of the orientation angle with increas- 
ing Win, seen in experiments for all samples is in qualita- 
tive agrement with the predictions of the DEMG model, 
resulting from the gradual unwinding and straightening 
out of the tube. The relative positions of the x versus 
Wi R curves from top to bottom in Fig. ^ for the three 
samples are the result of a competition between the ef- 
fects of Nf. and A them. Having a lower value of A (see. 
Fig. [l0| ) prompts ||2^ the oricntational angle for PS82 to 
fall below those for the other two samples, but as can 
be seen in Figs. |l^ and 20, the effect of N^, is stronger 
in this case too, since the curves with a lower value of 
(that for PS82, here) are predicted |2^ to show a 
lower degree of orientation at any fixed Wir in this flow 
regime. Similarly, PS2 exhibits a higher degree of orien- 
tation than for PS81 at a fixed Wir, since A is lower for 
PS2 in this intermediate regime of flow. As Wir — > 1, 
the tube is expected to straighten out to its full length 
and become orientated in the direction of the outflow axis 
[Fig. 1^ with the orientation angle reaching its asymptotic 
value, x', given in Eqn. (jl^). Using the asymptotic val- 
ues of A ~ 0.1 for PS2 and PS81 and A ~ 0.05 for PS82 
from Fig. |l0|, we obtain x' ^ 17.55° for PS2 and PS81 
and x' ~ 12.60° for PS82 which are in agrement with 
the results in Fig. |2^. Interestingly, the theoretical pre- 
dictions for the orientation angles for PS2 and PS81 in 
Fig. ^ show an "undershoot" before reaching x' ■ For 
the highest rates of deformation, Wir ^ 1.5, studied for 
PS82, the predicted orientation angle has not reached 
its asymptotic value. The experimental curves of x, as 
seen in Fig. qualitatively follow the behavior shown 
by the predicted curves except for that they do not show 
an undershoot behavior. In addition, x' ~ 17.5° and 
x' ^ 16° for PS2 and PS81, respectively and although 
the final value of x at the highest Wir studied for PS82 
is ~ 17.4°, it does not seem to have saturated; and for 
0.3 < WiR < 1.5 the x values for PS81 and PS82 overlap. 

The third flow region of highly non-linear viscoelastic- 
ity {WiR > 1) is predicted to show three distinct 
signatures. Firstly, the tube orientation should be com- 
plete as indicated by x = x'- We have discussed about 
the behavior of the predicted and experimental orien- 
tation angles for all polystyrene samples in the preced- 
ing paragraph. Secondly, the onset of chain-stretching is 
predicted to take place with a clear evidence of the di- 
mensionless birefringence exceeding its plateau value of 



unity as shown in Fig. This is evident in the experi- 
ments too, as displayed in Fig. |l^. Thirdly, the onset of 
chain-stretching is also marked by the collapse of all the 
birefringence curves to a single universal curve at the end 
of the plateau i.e., they should follow the scaling behavior 
An/ (CG'j^) = f{WiR) in this region, as is clearly shown 
by the predicted traces in Fig. |l^. This feature is well 
demonstrated in the experiments (Fig. |l^) by PS81 and 
PS2 but not by PS82, which fails to show a well-defined 
plateau too. At high rates of deformation {WiR^ > 1), the 
chain-stretching dynamics prompts the different birefrin- 
gence curves to reach different asymptotes proportional 
to ~ Tit/Ne- The theoretical curves in Fig. |l8| show these 
features: the increase in the birefringence with WiR in 
this regime is faster for PS82 (n* = 8420 and - 7) 
than for PS81 {nt = 8420 and 13), which is in turn 

faster than for PS2 (n* = 2890 and ~ 13). However, 
none of the An/(CG^) curves are predicted to reach 
their asymptotic values (~ 648 for PS81, - 1203 for PS82 
and ~ 222 for PS2) for the highest WiR studied, though 
the curves for PS81 and PS82 indeed show a tendency 
towards saturation in Fig. ^ At the highest rates of 
deformations, the experimental birefringence levels (Fig. 
p^ ) for all samples are significantly below their satura- 
tion values, and they do not show a tendency to saturate. 
Contrary to the prediction, the experimental An/{CG%) 
curve for PS81 falls below that for PS2. We note that the 
relative level of birefringence for the different samples as 



depicted by the individual curves in Figs. 17 and 18 are 
strongly dependent on the choice of the value of the relax- 
ation modulus, G% [Eqn. (^]. However, we believe that 
the qualitative as well as quantitative agreement between 
the experimental and predicted birefringence curves well 
justifies the choice of value for G^. 

Following Mead et al. [^, we define the generalized 
viscosity function (which is the conventional viscosity ap- 
propriately generalized for mixed shear and extensional 
flow) as follows: 



7(1 + A)- 



(19) 



This reduces to the usual deflnition of the shear viscosity 
and the planar extensional viscosity in the limits of A = 
and A = 1, respectively. The scaling factor (1 -f A) makes 
the viscosity independent of the flow-type A at low strain- 
rates. In the principal axes of the rate of deformation 
tensor, which are at 45° from the coordinate axes used 

in Figs. I and |, we get 77(7, A) = ^/f^TW- Using the 
stress-optical relationship [^ , which states that the bire- 
fringence tensor is proportional and coaxial to the stress 
tensor, one can relate the birefringence, which is the dif- 
ference of the principal values of the refractive- index ten- 



sor, n, in the plane of the flow (i.e.. An = n'^. 



the stress as Uxy 
written as 



2C 



sin(2x), so that Eqn. (19) can be 



vy 



), to 
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'7(7, •^) = 



An sm{2x) 
2C7(1 + A)' 



(20) 
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DEMG model predictions 




10"' 10"- 10"' 10° 

WiR 

FIG. 21. The predicted behavior of the steady-state gen- 
erahzed extensional viscosity scaled -with the zero-shear vis- 
cosity, as a function of Win for the polystyrene solutions, 
obtained by using the computed values of An and x from 
the DEMG model, the measured values of -y and A, the 
stress-optical relationship, and the definitions of rj and 770. 
The model parameters nt, N^, and tr are given in Table I. 

We present the predicted and the experimental steady- 
state viscosity function versus the measured Wiji for the 
three polystyrene solutions in Figs. and |2^, respec- 
tively. The predicted values of rj are calculated from 
Eqn. ( pO| ) by using the predicted birefringence and ex- 
tinction angle sho"wn in Figs. ^ and and by scaling 
the viscosity "with the zero shear viscosity. 



45 



(21) 



Similarly the experimental plots sho"wn in Fig. ^ are cal- 
culated "with the use of the measured An (Fig. 17) and 
X (Fig- ^) and then by normalizing by the experimental 
ryo given in Table |. In sharp contrast to the observa- 
tion that dilute solutions, (e.g., Boger fluids) sho"w 
strain-rate thickening at high strain-rates, for entangled 
solutions the viscosity function is predicted to sho"w dif- 
ferent behaviors in the three earlier defined flo"w regimes. 
In the first flo"w regime of linear viscoelasticity, the vis- 
cosity is predicted to be nearly constant and equal to 
its zero shear-rate value, r/o, (Fig. ^ ). This is approxi- 
mately follo"wed by the experimental curves (Fig. |2^ ). In 
the intermediate flo"w regime [bet"ween Wid > 0(1) and 
Wiu < 0(1)], where the stress is generated primarily 
via orientation of the tube segments, the model predicts 
a strong strain-rate thinning behavior "with the degree of 



thinning being an increasing function of and a de- 
creasing function of A, as sho"wn in Fig. ^ In the ex- 
periments, "we see that has a stronger effect than A 
on the degree of thinning. The predicted effect of A on 
the relative positions of 77/770 versus Win curves for PS81 
and PS2 from top to bottom is not supported in the ex- 
periments. 



Experimental results 




WiR 

FIG. 22. The steady-state generalized extensional viscos- 
ity scaled -with the zero-shear viscosity, versus Win for the 
polystyrene solutions, obtained by using the experimentally 
measured values of An, 7, A, rjo, the stress-optical rela- 
tionship, and the definition of rj. The data are plotted in the 
same scale as the corresponding DEMG model predictions in 
Fig. 21. 

The presence of an "undershoot" in the predicted ori- 
entation angle for PS2 (Fig. |2^) in this flo"w regime clearly 
manifests itself in the predicted viscosity plot (Fig. |2|) 
too. This does not seem to be the case "with PS82, be- 
cause for PS82 the small undershoot present in the pre- 
dicted orientation angle does not sho"w up in the pre- 
dicted viscosity. Although, the experimental viscosity 
function sho"ws a thinning behavior "with the degree of 
thinning being much "weaker in PS82 than the other t"wo 
solutions, the model predicts much stronger viscosity 
thinning behavior than observed in the experiments (ex- 
cept for PS81 "where they are similar). Also, the expected 
"undershoot" in the viscosity for PS2 is absent in the ex- 
periment. 

In the third regime of flo"w [Win — 0(1)] the model 
predicts a sharp upturn in the viscosity function as chain- 
stretching is predicted to take place. As the rate of de- 
formation is increased, the viscosity is also expected to 
increase and finally reach an asymptote for Win >> 1. 
Surprisingly, ho"wever, except for a small upturn sho"wn 
by the last three data points for PS2 in Fig. E2, these pre- 
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dieted features are completely absent in the experimental 
data. The cause for this lies in the fact that even though 
the experimental birefringence values of these samples 
(Fig. |l^ rise above the plateau value, indicating chain- 
stretching, the rise is much weaker than predicted (Fig. 
p^ ). Also, the degree of orientation of the tube segments 
in the experiments (Fig. |l^) fails to show the sharpness 
predicted by the model (Fig. EQ) in this flow regime. 
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FIG. 23. A comparison of the Weissenberg number de- 
pendence of the measured steady-state orientational angle, 
X (in degrees), with values predicted via the DEMG model 
using the measured values of both 7 and A for the PS81 so- 
lution. The model parameters were Ne = 13, rit — 8420, and 
TR = 2.25 s. 

The above two effects compete with each other to nul- 
lify any signature of chain-stretching that could have 
been otherwise seen in the generalized viscosity func- 
tion. We have noted before that the present model [ p3| , p4[ 
has improved the original DEMG model by using finitely 
extensible freely jointed chains instead of infinitely ex- 
tensible Gaussian chains. In the limit of the high value 
of Tit (or Muj), the original version of the model is re- 
trieved and at high deformation rates the chains stretch 
tremendously to a show a nearly singular viscosity or 
birefringence. This effect is apparent in the model pre- 
dictions shown in Figs. |2l| and 18, where both PS81 and 
PS82 (higher nt or M^) show a steeper increase in rj and 
An/ {CG'^) at Win ^ 1) compared to PS2, but is absent 
in the experiments (Figs. ^ and [l7|) . Unfortunately, the 
highest dimensionless rate of deformation or Win that 
we could reach in the experiments is only about 3 for 
PS81 and about 1 for the other two solutions. This defi- 
ciency is caused by the fact that the viscoelastic modifi- 
cation relative to the Newtonian flow was much stronger 
than we had initially expected. The measured values of 



7 and A were thus reduced compared to their Newto- 
nian values, so that the originally calculated theoretical 
(WiR)^^^ = (jNy/XNTR)^^^ = 0(10) corresponding 
to the maximum motor speed chosen for the experiments 
ultimately produced only (Win)^^^ = (-yVXTn)^^^ = 
0(1). The absence of chain-stretching signatures in the 
viscosity may also have to do with the choice of sample 
compositions which gives rise to low values for the num- 
ber of entanglements per chain. It would be interesting 
to carry out further experiments on samples with higher 
Nf, in this highly non-linear viscoelastic flow regime at 
high enough Win, particularly to look into the chain- 
stretching effects much more carefully. 
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FIG. 24. The experimental data for the steady-state, x (in 
degrees), versus Wijt for the sample PS82 compared with 
the corresponding DEMG predictions, obtained by using the 
measured values of both 7 and A, and the model parameters 
= 7,nt= 8420, and tr = 3.01 s. 

Let us now turn our attention into more quantitative 
comparisons between the steady-state results and the 
predictions of the DEMG model. We will first consider 
the case of the orientation angle % for three samples as 
shown in Figs. |2^, |2j, and Since there is no other 
parameter involved for scaling (e.g., in the case of the 
birefringence), direct comparisons between the experi- 
mental data and model predictions are possible in this 
case. Although the experimental data follows a similar 
qualitative decreasing trend in the orientation angle with 
increasing Wifj, as predicted, the quantitative mismatch 
with the prediction is quite obvious in these figures. In 
the limit Wiji 0, it is expected that x should be 45°, 
corresponding to alignment of the principal optical axis 
with the principal axis of the rate of strain tensor. The 
experimental accuracy of both An and x s-t the lowest 
rates of deformation is limited by the sensitivity of the 
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apparatus in determining a minimum birefringence and 
its associated orientation. Also, the contribution from 
the residual glass birefringence is important at these low 
levels of polymer anisotropy. At low and intermediate 
Win the DEMG model underpredicts the orientation but 
at higher Win it overpredicts the same, and there is a 
crossing point between these two behaviors. By com- 
paring the results on the orientation angle with those 
obtained from the DEMG model numerical simulations 
with identical parameters but with a constant flow-type 
parameter A = Xn = 0.1501, we have confirmed that 
the opposite curvatures seen in the predicted curves be- 
low and above the crossing point depends on the way A 
changes with Win. The fact that the DEMG model over- 
predicts the tendency of the flow to orient polymer chains 
towards the outflow axis at intermediate and high strain- 
rates has been observed in earlier experiments |l^ , [l^ and 
it was speculated that "tube-dilation" , which is not in- 
corporated in the DEMG model, may be responsible for 
such an effect. In simple terms, this means that the ef- 
fective dilation of the tube radius is stronger for large 
rates of deformation, leading to a reduction of and 
hence the lower orientation angle seen in the experiments. 
Weather or not tube-dilation is present in the system, for 
Win > 0(1), we expect that the polymer molecules will 
become oriented along the outflow axis of the flow-field. 
As we have noted earlier, and as can be seen clearly from 
these figures, the asymptotic value of x is approximately 
reached in experiments for PS81 and PS2, but not for 
PS82. 
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FIG. 25. A comparison of the Weissenberg number depen- 
dence of the measured steady-state orientational angle, x (in 
degrees), with values predicted via the DEMG model using 
the measured values of both 7 and A for the PS2 solution. The 
model parameters were A^e = 13, iit ~ 2890, and tr = 0.56 s. 
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FIG. 26. A comparison of the measured steady-state scaled 
birefringence, An/{CG%) at several different Win, with the 
values predicted via the DEMG model obtained by using the 
measured 7 and A for the PS81 solution. The model param- 
eters used were — 13, rit — 8420, and tr = 2.25 s. 
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FIG. 27. The experimental data for the steady-state 
An/{CG%) versus Win for the sample PS82 compared with 
the corresponding DEMG predictions, obtained by using the 
measured values of both 7 and A, and the model parameters 
= 7, nt= 8420, and tr = 3.01 s. 

The qualitative as well as quantitative match between 
the experimental and predicted birefringence curves for 
these entangled samples, shown in Figs. |2^, and 28, 
is quite satisfactory. This firstly points to the fact that 
the choice of the plateau modulus, as described earlier 
in this section, works very well for these entangled sys- 
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terns. Given the proper values of G^, we can clearly 
see that DEMG model predictions for birefringence are 
fairly close to the experimental measurements, at least 
for low and intermediate Wi^. We note that the plateau 
expected in the transition region between the dynam- 
ics dominated by segmental orientation and segmental 
stretch, is quite narrow, because of the small separation 
between Td and {tii/tr = BNg = 21 for PS82, and 
39 for PS81 & PS2). The plateau is smeared out in the 
experimental curves, since tube-dilation reduces and 
hence the separation 3A^e between these two time scales. 
This supports the tube-dilation idea. In the nonlinear 
viscoelastic regime, the chain-stretching predicted by the 
model is much stronger than what has been seen exper- 
imentally. In the limit Win — > oo, the model predicts a 
saturation of birefringence ^ nt/Ne corresponding to a 
maximum chain-extension ratio ~ y^nt/Ne- Comparing 
the maximum birefringence shown in Figs. |2^, and|2^ 
with the above numbers, we see that the maximum chain- 
extension predicted by DEMG model at the highest Wir 
studied are 87.8%, 63.2% and 28.5% for PS81, PS82 and 
PS2, respectively but the experimentally observed max- 
imum chain-extension for these samples is only 15%. 
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FIG. 28. A comparison of the Weissenberg number depen- 
dence of the measured steady-state birefringence, An/ {CG%), 
with the corresponding DEMG model predictions obtained by 
using the measured values of both 7 and A for the PS2 solu- 
tion. The model parameters were A'^^ = 13, nt = 2890, and 
TR = 0.56 s. 

The quantitative agreement between the experimen- 
tal and the predicted rj/riQ for small rates of deforma- 
tion in Figs. and ^ directly points to a proper 
choice of the calculated and the experimental values of 
riQ. In the intermediate flow regime the DEMG model 
predicts a strong strain-rate thinning behavior for the 



generalized extensional viscosity for entangled solutions. 
This also distinguishes the entangled solution from di- 
lute ones, for which a monotonic increase in rj with the 
increase of 7 is expected. In the experiments, we see 
that the solutions with iVg ^ 13 show much stronger 
thinning than the solutions with a lower value of TVe, 
consistent with the prediction. It has been shown in ear- 
lier experiments that the original DE model overpredicts 
the strength of shear-thinning in simple shear flows jl3j 
and the DEMG model also shows a similar behavior in 
extension-dominated flows p^ . Similarly, we note that 
the excessive shear-thinning predicted by DEMG model 
for PS82 and PS2 (shown in Figs. |3^ and respectively) 
is due to the fact that the predicted orientation angle and 
the birefringence are too small in this intermediate range 
of WiR (see. Figs. H, ||, ^ and H). The tendency of 
DEMG model to over-orient the polymer chain away from 
the principle axis of the rate of strain tensor, is specu- 
lated |p^ , p^ to be due to the fact that the model does not 
incorporate the tube-dilation effect which may be present 
in the real systems. For PS81, as can be seen from Fig. 

, the experimental curve for the viscosity closely fol- 
lows the predictions of the DEMG model in the low and 
intermediate range of WiR primarily due to the fact that, 
in this region. An is too small but x is too large so as to 
compensate each other. It can be seen from these figures 
that as WiR approaches 0(1), the model predicts an up- 
turn of the viscosity followed by a sharp increase owing to 
the prediction of a strong chain-stretching phenomenon. 
Apart from a small upturn seen in Fig. |3l|, the model 
fails to account for the observed behavior of the gener- 
alized extensional viscosity in this non-linear viscoelastic 
regime. 

1. Cox-Merz superposition 

The empirical Cox-Merz relation states that the steady 
shear viscosity is equal to the modulus of the complex dy- 
namic viscosity evaluated at the the angular frequency 
equal to the shear-rate. This is observed to be valid 
mostly in the linear viscoelastic regime [Q, and also in 
the case of a simple shear flow (A = 0), primarily because 
chain-stretching is not significant in either of these cases. 
Thus, for extensional flows (A > 0) departure from Cox- 
Merz "rule" constitutes another means to study chain- 
stretching effects present in a given flow system. Fol- 
lowing a suggestion of Mead et al. in Ref. |24| , we have 
compared the measured generalized extensional viscos- 
ity, T], for the steady-state flow generated in a two-roll 
mill as a function of the measured velocity-gradient (or 
"shear-rate") 7 for the entangled solutions with the cor- 
responding dynamic linear viscoelastic measurement, i.e., 
|7/*(ti;)|, in the form of "Cox-Merz plots" shown in Figs. 
|29| , pO| , and ^ Both parameters are scaled with the zero- 
shear viscosity tjo (see. Table |). The match between the 
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two sets of data at low deformation rates for all samples 
and that the values for these parameters are in the same 
ballpark provides an increased confidence on the choice 
of the corresponding ryo values. At low and intermediate 
range of flow deformations, the values of these two sets of 
parameters arc in a closer agreement for the samples with 
TVe 13 than that for PS82 (TVg ~ 7). Apart from this, 
the overall dissimilarities between the deformation rate 
dependence of these two sets of results is obvious from 
these plots. At high deformation rates, chain-stretching 
effects become important, and, though it is not as strong 
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as demanded by the DEMG model seen in Figs, 
and ^ in the previous section, its effect is strong enough 
to show an increased departure of i] from the Cox-Merz 
superposition. This departure is quite expected, as we 
have noted above, since the chain-stretching effects are 
not accounted for in this empirical rule. The departure 
seen in the intermediate range of deformation rates is 
because of the fact that even in this regime, the flow ex- 
perienced by the polymers in the stagnation region of a 
two-roll mill is very different from a shear- flow. In short, 
these figures provide a quantitative picture of the relative 
importance of chain-stretching to orientation effects. 
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FIG. 29. The measured steady-state generalized exten- 
sional viscosity for the PS81 solution scaled with the 
zero-shear viscosity at several different Win is compared with 
that computed from the DEMG model. The experimental re- 
sults were obtained from the measured values of An, x, 7, A, 
770 , and using the stress-optical relationship and the definition 
of 77. The theoretical values were obtained from the computed 
values of An, x and 770 using the DEMG model with the mea- 
sured values of 7 and A as input. The model parameters were 
Ne = 13, nt = 8420, and tr = 2.25 s. 
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FIG. 30. A comparison of 
scaled generalized extensional 
the PS82 fluid (obtained from 
7, A, 770, and with the use of 
and the definition of 77), with 
from the computed An, x a^ud 
(iVe = 7, nt = 8420, and tr = 
and A as input. 
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the experimental steady-state 
viscosity 77/770 versus Wir for 
the measured values of An, x, 
the stress-optical relationship 
the predicted values obtained 
770, by using the DEMG model 
3.01 s.) with the measured 7 



IV. SUMMARY AND CONCLUSIONS 

The recently developed dynamic light scattering 
section of the two-color flow birefringence experimental 
setup is used to characterize the steady-state flow-fields 
for a Newtonian fluid as well as for viscoelastic, entan- 
gled polymeric fluids in a two-roll mill, by measuring the 
velocity-gradient and the flow-type parameter for these 
fluids. Assuming Newtonian, creeping flow symmetry of 
the flow-field, the flow parameters for a Newtonian fluid, 
obtained directly from the decay rate of the autocorrela- 
tion functions measured at the stagnation region of the 
flow, are shown to be very well described by the Frazer's 
creeping flow solution . 

The measured flow parameters for the polymeric fluids 
show clear departures from their corresponding Newto- 
nian values, providing direct evidence of flow modifica- 
tion due to the conformational dynamics of the polymer 
molecules. Our results suggest that this departure may 
have a stronger dependence on the entanglement density 
of the polymer chains than on the molecular weight and 
concentration. 
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FIG. 31. A comparison of the steady-state generalized ex- 
tensional viscosity for the PS2 solution scaled with zero-shear 
viscosity versus Wir, and that computed from the DEMG 
model. The experimental results were obtained from the mea- 
sured values of An, x, 7, A, rjo, and by using the stress-optical 
relationship and the definition of r]. The theoretical values 
wore obtained from the computed values of An, x and rjo, by 
using the DEMG model with the measured 7 and A as input. 
The model parameters used were Ne = 13, nt = 2890, and 
tr = 0.56 s. 
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FIG. 32. The modulus of the complex viscosity, |77*|, versus 

the angular frequency, uj (in rad/s), for the polystyrene solu- 
tion PS81 is compared with the steady-state generalized ex- 
tensiorial viscosity, rj, versus the measured velocity-gradient, 
7 (in s~^): the "Cox-Merz plot". Both results axe scaled by 



the measured zero-shear viscosity. 
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FIG. 33. A comparison of \ri*\ and 77 for the sample PS82 
plotted against the radial frequency, cu (in rad/s), and the 
measured velocity-gradient, 7 (in s~^), respectively: the 
"Cox-Merz plot". Both results are scaled by the measured 
zero-shear viscosity. 
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FIG. 34. The "Cox-Merz plot" for the polystyrene PS2 
fluid, i.e., a combined plot of the scaled modulus of the com- 
plex viscosity, |t;*1/?7o, versus the angular frequency, lo (in 
rad/s), and the scaled steady-state generalized extensional 
viscosity, r?/??o, versus the measured velocity-gradient, 7 (in 
s-i). 

Within the hmit of experimental error, the flow-field 
at the stagnation region of the t"wo-roller is verified to 
retain its symmetry, for all rates of deformation studied 
with the Newtonian fiuid and for the strain-rates 7jv\/Ajv 
< 0.5 s~^ {Win < 1) with the polystyrene samples. The 
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flow-type parameter for the polystyrene fluids, extracted 
from the steady-state DLS experiments, was found to 
exceed its Newtonian value at very low rates of deforma- 
tion. This result was verified to be consistent from three 
different experimental means. 

The dynamics of the entangled polymers under study 
was coupled with the changes in the flow-field. As we 
have noted above, the DLS experiments has probed the 
effect of the viscoelasticity on the velocity field compared 
to its Newtonian form. On the other hand, the poly- 
mer response induced by the flow-fields generated by a 
two-roll mill is studied using two-color flow birefringence 
experiments, which were in turn compared with the bire- 
fringence predictions from the DEMG model using the 
measured flow data as input to the model. 

Similar to the model predictions, the experimental 
birefringence results clearly illustrate the existence of 
three steady-state flow regimes. In the first two regimes 
of low and moderate rates of deformation, where the dy- 
namics is dominated by chain segment re-orientation and 
reptational diffusion, the model predictions are qualita- 
tively, and some times, quantitatively, reproduced in the 
experiments. However, the model is clearly inadequate in 
describing the polymer dynamics at sufficiently high rates 
of deformation, where chain-stretching is important. Our 
data do show signatures of chain-stretching: the experi- 
mental birefringence exceeds the plateau value, the Weis- 
senberg number scaling behavior shown at the onset of 
chain-stretching by at least two of the three polymer sam- 
ples studied, and the departure of the generalized exten- 
sional viscosity data from the empirical "Cox-Merz super- 
position" . The comparisons of the experimental birefrin- 
gence as well as the generalized extensional viscosity data 
unambiguously indicates that the DEMG model overpre- 
dicts chain-stretching in these entangled samples. The 
relative values of birefringence, orientation angle, and 
viscosity for different entangled polymeric samples are 
shown to be determined by the competition between the 
effects of the corresponding flow- type parameter, molec- 
ular weight (via nt) and the number entanglements per 
chain. Our experiments have demonstrated that the ef- 
fect of Ne on these parameters is much stronger than that 
of the AI^ and A. We note that the model fails in describ- 
ing the "smearing" of the expected plateau region in the 
birefringence curves. In particular, it overpredicts the 
tendency of the flow to rotate the polymer chains toward 
the outflow axis, at high Weissenberg numbers, thereby 
predicting excessive thinning of the generalized viscosity 
in this regime of flow for at least two of the polystyrene 
samples of our study. It was suggested in an earlier pa- 
per [|l2| that one of the prime causes for these deflciencies 
may be the fact that a real system experiences a confor- 
mation dependent decrease in the entanglement density 
(or a dilation in the tube diameter) as well as convective 
constraint-release, causing a decrease in the time scale 
for reptation with increasing flow strength that is not 



included in the present form of the DEMG model. This 
calls for further efforts to improve reptation based consti- 
tutive models to obtain a better match of its predictions 
to experiments on entangled polymeric systems. 
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TABLE IIL The values of the parameter h{(f)) for three 
orientations 4> of the two-roll mill and that for the special 
case of a symmetric flow (e = 0). 



<p 


General 


Symmetric flow 


0° 




7^ 


90° 




f A^ 


45° 


e2 + i7'(l + A2) + e7(l-A) 





TABLE I. The cliaraetoristic: paranictcrs for the throe 
polystyroiu' samples. 

Sample c 770 A^^ tr n,t K n 

(xlO^) (^) (P) (s) [P(^)"] 

PS81 8.42 0.0396 7500 ~ 13 2.25 8420 343.72 -0.80 

PS82 8.42 0.0262 2700 ~ 7 3.01 8420 190.11 -0.70 

PS2 2.89 0.0867 19000 - 13 0.56 2890 2703.96 -0.83 



TABLE IL The characteristic parameters for the two-roll 
mill with radius R, height H, and gap G. 

Roller R (cm) G (cm) # (i^)exp {Xn\^ (Aiv)exp 

G 1.278 0.844 3.1 2.7 2.8 1.501 0.1508 
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